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SUMMAEJY 


On the hasis of a recently developed, theory for finite swept- 
hack wings at supersonic speeds, calculations of the supersonic 
wave drag at zero lift Trare made f cr a series of ^flngs having thin 
symmetrical biconvex sections with untaperod plan foamis and various 
angles of sweephack and aspect ratios. Tne results are presented 
in a unified, form so that a slng.le chart permits the direct deter- 
mination of the \rave drag for this family of airfoils for an exten- 
sive range of aspect ratio and sweephack aaagle for stream Mach 
nvimhers up to a value corresponding to that at which the Mach line 
coincides with the wing leading edge. 

The calculations showed that in general the wave-drag coeffi- 
cient decreased with Increasing sweephack. At Mach numbers for 
which the Mach lines are appreciably ahead of the wing leading 
edge, the wave-drag coefficient decreased to an important extent 
with increases in aspect ratio or slendomess ratio. At Mach 
numbers for which the Mach lines appx’oach the ■vdng leading edge 
(Mach numbers approaching a value equal to the secant of the angle 
of sweephack), Ihe wave-drag coefficient decreased with reductions 
in aspect ratio or slenderness ratio. 


INTOODUCTION 


Recent developments in airfoil theory for supersonic speeds 
(references 1 and 2) indicate pronounced effects of sweephack and 
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aspect ratio on the dr'ag. In reference 1, a theory vae developed ^ 

for calculating the pressure distribution at supersonic speeds and 

zero lift for swept-hack wings of arbitrary linear taper and aspect 

ratio having thin symmetrical airfoil sections of sharp leading 

edges . 

In the present paper, the method of reference 1 is applied to 
calculate the supersonic wave drag for a series of wings having 
thin symmetrical biconvex sections at zero lift with untapered 
plan forms and various angles of sweepback and aspect ratios. The 
term "biconvex profile" as used herein refers to an airfoil secticaa 
composed of two parabolic arcs. In each case, the wing is con- 
siderod to be cut off in a direction parallel to- the direction of 
fli^t. In the calculations the Mach number is varied from 1.0 
to a value corresponding to that at which the Mach line coincides 
with the wing leading edge. The results of the calculations are 
presented in a unified fonn -v-iilch permits the direct determination 
from a single chart of the wave drag for this family of wings for 
an extensive range of sweepback angle and aspect ratio for Mach 
numbers from 1.0 to the value coxrespondlng to that at which the 
Mach line coincides with the wing leading edge, or eq.ual to the 
secant of the angle of- sweepback. Although the calculations have 
been made for the biconvex profile, the data may be applied to 
indicate corresponding results for profiles approximately similar » 

to the biconvex. 


SYMBOLS 


X, y, z 


dz/dx 

c 

t/c 

A 


coordinates of mutually perpendicular system of axis 
in wing 

slope of airfoil surface 

chord of airfoil section, measured in flight direction 
thickness ratio of section, measured in fli^t direction 
angle of sweep, degrees 


m = cot A 


r 



NACA TS No. 1319 


3 


h vlng semlspan measured along y-axis, seamLciiords except 

In appendix A 

K parameter indicating span-wise position equal to j/m., 

semi chords 


A 

l/t 

P 

p/g. 

T 

u 


u 

V 

I 

M 


aspect ratio 

slendemesa ratio, ratio of -idLng semispan measured along 
leading edge to maximum thickness of center section 

disturhance pressiure 

pressure coefficient, ratio of disturhance pressure to 
dynamic pressure in free stream 

velocity in flight direction 

x-cong)onent of disturbance velocity, positive in flight 
direction 

u caused by source line -with reversal in sign of m 

z-coarponent of disturbance velocity 

disturbance-velocity potent 5 .al 

source factor required to maintain a given vedge angle 

Mach number 


P = - 1 


d 

°d 

Ac, 


o& 


coordinate measijred along y-axis -which is shifted to tip 
section, semichords 

coordinate measured along y-axis -vdiich is shifted to 
opposite tip section, semichords 

•wave drag at section 

-wave-drag coefficient at section -id-thout tip effect 

■wave-drag coefficient at section including tip effect 

increment in section -wave-drag coefficient ca'used by -wing 
tips 
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Ac^ Increment in section wave-drag coefficient caused by wing 

^ tip located on same half of wing as section 

Ac^ increment in wave-drag coefficient at section on one wing 

■ panel caused by tip of opposite wing panel 

Cp wave -drag coefficient for wing without tip effect 

00 

Cp wave-drag coefficient for wing including tip effect 

ACj) increment in wave-drag coefficient caused by tips, com- 

plete wing 

^p increment in wave-drag coefficient on one wing panel caused 

by adjacent wing tip, complete wing 

ACp increment in wave-drag coefficient on one wing panel caused 

by tip of opposite wing panel, complete wing 

I, Ti auxiliary variables which replace x and y, respectively, 
used to indicate origin of source line 

Primed values of A, y, y^, y^, h, m, dz/dx, and z indicate 

transformation involving multiplication by factor p. 

Subscripts; 

1, 2 wings related according to transformation which malces yP 
and mp equal respectively for two wings 

Subscript notations for u and u indicate the origin of source 
line in terms of coordinates x and y, respectively. 


ANALYSIS 


Basic data .- The present analysis is based on thin-airfoil 
theory for small, pressure disturbances and a constant velocity of 
sound throughout the fluid. The axes used are the mutually per- 
pendicular X, y, z system in which the x-axls is taken in the 
direction of fll^t positive to the rear, the y-axis is along the 
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span positive to -Uie rif^t, and the z-axts is positive up^wards. 
Figvire 1 shovs the syaibols used to desiGnate the -VTing-plan-fona 
parameters. The present analysis is made for untapered -wings of 
■biconvex profiJ.e at zero lift and is limited to a Mach number 
range from 3..0 to -tlie value corresponding to that at -vdiich the Mach 
line coincides -with the wing leading edge^ that is, to a value 
equal to the secant of -fche angle of sveepbaclc. 

Theory ." If p is -the distfirbance pressure conrputed for one 
surface of the airfoil section, the -wave drag for -the section is 


d = 2 




( 1 ) 


and the section -wave-drag coefficient 


°d 



dx 

q. dx 


( 2 ) 


"where dz/dx is the slope of the surface of "the airfoil at "the 
point X. For the symmetrical biconvex profiles (ccanposed of 
parabolic arcs) 


dz 

dx 



-where the thickness ratio t/c may be considered tho sole shape 
parameter. From thin-airfoil theory. 


p 2 u 

q. V ox V 
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■where u Is -the distrrhance-velocity ccaaponent in the x-diroction^ 
taken positive in the fli^t direction. For a given S'wept-’baclc 

■srtng with A a cot the di’ag coefficion’fc-in equation (2) at 

the span'wise station 7 may no'vr he ■va'itten as 


Cfl(y) - 


„2 


n|fc 


y/m 




(3) 


The desired integrand u in equation ( 3 ) may he determined 
hy the procedure given in reference 1. On the oasis of the 
linearized theory, reference 1 derives a solution representing an 
oblique ( swept -hack) sounce line making the angle of s'weophack A 
with the y-axis. This soluticsn utilized for the pressure field 
or for -bhe disturhance velocity is the real part of 


Uq,q « I cosh 


1 X - 
P |y •• mx { 


(4) 


where the subscript notation indicates that the source line starts 
at -bhe origin of coordinates (x = 0, y = O) • Equation (^v) is 
shown in reference 1 to satisfy the boundary condition for a thin 


oblique ■wedge mak'ing tho half -angle 


coordinate syst em of re forenco 1 

/dz Y _ _ 5, Vi '- 

\dx/ V " V m' 


(y' 


( dz Y 

/ 

yP^ 


idx - 


in the transformed 
z * = zp) where 


I. In order" to "obtain an equal 'vrodgo angle 


in the physical coordinate system, the following relations between 
the transformed coordinatos of reference 1 and "^e physical coor- 
dinates are used 


m' = mP 


■w' = 



w = 


Bz 


r 
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■where w* and w refer, respectively, to •the vertical velocities 
in -the •transformed and physical coordinate systeaas. Tlius 


5z * 5z 


33 


w*p 


The half-angle of the wedge is then determined as 


^ w It \jl - 

dz Y V m 


or the req.uired source factor in order to maintain the desired 
wedge angle is 


T m dz 

yi - 


.(5) 


By superposi'triLon of solutions of "the wedge type, swept -hack 
wings of desired profile shape and plan form, can he huilt up 
(reference 1) . In order to satisfy the boundary conditions over 
■fche surface of a biconvex wing, semi -infinite source lines of 
eq.ual s-trengfch are placed along the leading and trailing edges 
beginning at the center secticm, in conjunction with a constant 
dis-fcribution of sink lines along the chord. At the tip, ifiiere the 
wing is cut off, reversed semi-infinite source and sink lines are 
distributed' so as to cancel exactly the effect of those originating 
at the center secti<m in the entire region of space outboard of the 
tip. Tn "bhe present analysis, the tip is assumed to be cut off in 
the direction of fli^t. The term "tip effect" refers to the 
effect of this •wing cut-off. The form of the integrand u for 
eq,uation ( 3 ) is given in appendix A. 

In calculating the wave drag over the wing, the disturbances 
due to the elementanry sources and sinks are evidently limited to 
the regions enclosed by their Mach cones. Figure 2 shows the 
typical Mach lines originating from the source lines at the leading 
and trailing edges of the center and tip sections over a wing of 
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^5 sweepback. Figure 2(a) shovs the Mach lines for the infinitely 
long "wing^ and figure 2(b) includes the Mach lines starting from 
the tip section. In each case the disturbance over the ving caused 
by the leading- and trailing- edge source lines is limited to the 
region of the wing behind the corresponding Mach line. The regions 
affected by each of the Mach lines are indicated in figure 2(b) as 
regions I to IV . Region I represents the part of the wing affected 
by the sink line starting from the leading edge of the adjacent 
tipj region II represents the wing area affected by the sink lino 
starting from the opposite tip. Region III is influenced by the 
source line starting at the leading edge of the center section and 
includes the entire wing; region IV is Influenced by the source 
line starting at the trailing edge of the center section. 

The resultant velocity u at a point on the wing is made of 
the con^onent velocities caused by each of those source and sink 
linos where the influence of oach ccai5)onent is restricted to the 
ireglon behind its Mach lino. The drag coefficient Ca is there- 

fora obtained by evaluating in oq.uation (3) tho integrand Uq q 
over tho entire section (rogion III), tho integrand u^^q over 

part of tho section included in region IV, and the integrands for 
the u-oempononts caused by tho sink distributions along tho pro- 
file (fig. 2(b) and appendix A, equation (AS)). Tho drag coeffi- 
cients Aca and Aca _ are obtained simll^ly by integrating 

U.J 

along the section in tho rogions I and II, respectively, in addi- 
tion to the integrations for tho u-compraionts caused by the source 
distributions along tho profile (appendix A) . Tho limits of inte- 
grations for X along the chord and for y along the span, which 
represent the boundaries for tho rogions of influence for tho 
individual u-functions required fer a biconvex profile, are given 
in tho table in appendix A. 

Formulas for section wavo-drag coefficients .- The formulas 
for tho section drag coofflcionts obtained by intogration of the 
u-functions and by use of equation (3) ore presented in appendix B. 
Those foimulas give ospressione for the drag coofflciont without 
tho tip offect Cfl and also the expressions for tho increments 

in caused by tho tip effect Ac^^. 


Wavo-drag coefficients for complote wing .- In tho present 
investigation the section drag coefficients expressed by tho 
equations in appendix B woro Intogratod graphically to obtain tho 
results for tho wing-drag coof f iclonts . Subsequently, however. 
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analytical e3i)reasions for the IntegratlonB wore obtained. These 
formulas for the wing-drag coefficients are presented in appendix B. 

Drag coefficient of swept-back wing at Mach number of 1.0 . - 
The solution of the equations for Oj, given in appendix B shows 

tMt, for a symmetrical imtapered finite swept-back id.ng at Mach 
number of 1.0 and zero lift, positive and negative infinite values 
for are obtained at various sectioixs of the wing. The inte- 

gration over the wing of the limiting values for these infinite 
torms, howovor, gives zero. Althou^ sms sections of the wing 
have infinitely positive or negative drag, the total drag coeffi- 
cient over the wing results in a finite value. The prodiction of 
infinite vaJLuea of drag at certain sections of the wing clearly 
violates at theso sections tho assumption of small dlsturbancos 
from which the linearized theory is derived. The calculated values 
for the total drag coefficient at Mach number 1.0 are therefore 
guestionable. The formulas fcr the total drag at Mach number 1.0 
are presented in appendix B. 

Conversion of drag solution to series of related wings .- An 
examination of equations (Bl) , (Bj) , and (B5) indicates that the 
drag solution obtained for ono value of m and M may bo applied 
directly to obtain the drag for a whole sex*ies of wings in which 
each wing is at a certain appropriate Mach nuiiiber. (Eq.uation (Bl) 
is foimed by adding expression (Bib) to the right-hand side of 
egiiation (Bla).) For example, eciuation (Bl) may be expressed in 
tho following form; 



vdiore 



refers to tho variable terms and where mj3 = m'. 


and - = E. If tho subscript 1 refers to a -wing at Mach number 
m 

corresponding to and tho subscript S to any othei’ wing at the 

Mach number corresponding to Pg, then tho drag coefficients for 
the two wings may bo obtained from og.ua tion (6) as foU-ows; 



8 /t \ ^ %i^l 
\ ° / 2 _ Pi 


F 




(7a) 
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8 /t\ ^ J 72^2 

2 I moPp^ 


ni23. 


(713) 


Eqimtion (7) shovs tliat if = 721^2 “l^^l “ “2^2 


(■fe/c)^^ Pg ('l3/c)i^ cot A^ 

\ ^2 (t/c)g^ ^1 (t/c)g2 cot A2 


( 8 ) 


where c^ and c^ refac to tho spamd.ee positions and 


yiPi 

3o 


t respectively. In a similar manner, it can "be shown that if 


two wings are related according to y3_Pi = “l*^l " “2^£-’ 

the section drag coefficients obtained for wings 1 and 2 from equa- 
tions (B3) and (B5) aroL-in the same ratio as that eaproesod in 
equation (8). Equation (8), therefore, may he generalized to apply 
to the total drag coefficient at the section or, 


'dx = °d2 


(t/c)i^ cot Ax 


(t/c)g'^ cot Ag 


(9) 


where 

yiPi 
P2 ' 


and refer to the spanwlso positions and 

respectively. 


The wing-drag coefficients for wings 1 and 2 are given, 
respectively, by 





dyi 


(10) 


r 
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°d2 ^2 


( 11 ) 


By substituting for the integrand 


cdi 


in equation (lO) the rela- 


tionship expressed in eq.uation (9) and hy substituting 
for dy^_, equation (lO) may be written as: 


aygp2 

Pi 



(t/c)j_^ cot Ai 
h 2 (t/c) 2 ~ cot Aa yo 






or 


cot 

9d, = 9d ^ 

it/c)i£' cot Ag 


(It) 


Equation (12) permits a rapid determination of the drag coeffi- 
cient for winf^s of arbitrary sweepback, aspect ratio, and thickness 
ratio (within limitations of airfoil theory) from data obtained for 
one Bwept-back -v/ing for the appropriate range of aspect ratio and 
Mach number. For this purpose, use of a wing of 45° sweepback as 
the reference wing is most convenient. If the subscript 2 ie used 
to refer to the parameters for the wing of 45° sweepback and the 
subscript 1 is dropped, equation (12) becomes 


Cd = 


Cj)2 (V°)^ CO'*' ■A- 
(t/c)^® 


( 13 ) 


where and Cp 


2 


refer to wings idiose aspect ratios and Mach 


numbers are related by the following equations: 
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A2 a A tan 

A 

(14) 

Pg = p cot 

A 

( 15 ) 


The foregoing analysis shova that the results obtained for a 
vlng of I|-5° sveophack and a giT<aa aspect ratio Ag can bo trans- 
formed by moans of equation (13) to all vings for vdilch tho aspoct- 
ratio parameter A tan A = Ag and tho Ifech number parameter 

p cot A = ^2’ 

The grouping of the parameters as indicated in the foregoing 
axialysis permits tho use of a single genorallzod chart for pre- 
senting tho drag results. This chart is discussod in the section 
entitled 'iResults and. Discussion. " 

Prandtl-Glauert rixlo mod.ified for supersonic flo-ws .- In con- 
sidering the linearized problem of supersonic flov past a wing; it 
is often convenient to refer tho supersonic results for a given 
wing to a transformed wing at a rofeyence Mach n.umbor of M^ a 2, 

If this transformation is used a rule rosombling tlio Prandtl-Glauort 
rule (reference 3) for the aubsonio case, whervp M * 0 is 
the roforonce liach number, may be obtained. This rulo may be stated 
as follows* 


The streamline field of tho supersonic flow for a given body 
at a stream Mach number M may be calculated bj- multip.lying tho 
given y- and z-dimonslons, incl\iding those for tho Mach linos, by 

tho factor \/m^ - 1 and then by calculating the flow about the 


resulting transformed body at tho Mach number \fp.. Tho proBsm’e p 
and velocity incromonts u for tho given body ct tho Mach numbor M 
can then be obtained by multlplyiiog tho colculatod pressuro p and. 
velocity incromonts u at corresponding points of tho transformed 


body by tho factor 



It is interesting to note that the derivation of formulas (13) 
to (15) as given in this papor corresponds to utilizing tho solution 
for a transformed wing for tho •vdxolo family of wings rolatod to 
this transformation and tlien applying the aforomontionod modified 
Prandtl-Glauort rulo. 


r 
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Yarlatlon of section drag coefficient along span .- Figure 2, 
vhich was introduced previously to illustrate the system of Mach 
lines, also shows the variation of section drag coefficient Cjj_ 

along the span. The data, are presented for a wing of 45° sweephaclc 
and thickness ratio of 0.10. Figure 2(a) gives the results for a 
wing of Infinite aspect ratio and figtare 2(h) gives the results for 
a wing of finite span. 

In figure 2(a), the data are shewn for Mach nvmibera of 1.100, 
1 . 343 , and 1.4li(. The 3.owest Mach number 1.100 represents a case 
in Tdiich the wing leading edge diverges r-apidly from the Mach line 
(upper part of fig. 2(a)). In this case, tho section drag coeffi- 
cient c^ has a inaxiiirun value of 0.0542 at the center section, 
then drops sharply to zero at a distance of 1.13 chords from tho 
center line; Beyond this point, the wing shows a negative drag, 
idilch approaches asyuqptotically the suhsonic value of zero at an 
infinite distance from the wing center. This tjTpe of wave-drag 
distrihution is similar to that indicated in figure 11 of refer- 
ence 1 for a wing of 60° sweophack at a Mach nvanher of 1.4. 

For tho higher Mach numher 1.343, tho spanifiso variation of cg_ 

is markedly flatter. Unlike tho preceding caso, the drag coeffi- 
cient does not have its maximum value at tho center section "but, 
at first, increases in tha outboard direction, then roaches a peak 
and falls to zero at a distance from tho center of 6.6 chords 
(not shown in fig, 2(a)). 

At tho highest Mach number 1.4l4, the Mach line becomes coin- 
cident with the wing loading edge. In this case, tho wing gives 
a very hi^i drag and tho section drag coofficiont incroosoa in tho 
outboard direction, approaching Infinity at an infinite distance 
from the wing center. 

Figure 2(b) illustrates tho condition at which the aspect ratio 

is less than i/\^ - 1. Tho calcTilatod case shown is for an aspect 
ratio of 1.86 and Mach number of 1.10. In this caso, tho two wing 
tips cause incremonte in section drag coofficiont on each half of 
the wing, namely, and ^djj* affoct ^dj ^ 

given distance frm the tip is independent of tho aspect ratio. Tho 

tip effect AC(3 , however, is a function of tlie asnoct ratio. 

H 
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Effect of vlng tip on tdng-drag coefficient .- Eigure 3 Bho"vre 
the typical variatioa -with aspect ratio of the increment in Cjj 
duo to the tip ACJjj. The data are shown for a wing of 45° sweep- 

hack and a thickness ratio of 0.10 for tiach numhors of 1.100, 1*343, 
and 1.4l4. The present analysis for the untaperod wings indicated 
that if the aspect ratio is equal to or greater than 2m/(mp + 1), 
the integrated value of ^dj over the %d.ns is zero. On this basis, 

if the aspect ratio of the wing is greater than i/\/m^ - 1, the 
total increment in drag contributed by tip is zoro. As tho aspect 
ratio is reduced, the tip effect wiiich occurs when tho 

aspect ratio is smaller than - 1, however, loads to an 

increase in Op. The tip effect ZsCp then roaches a peak value at 

a certain aspect ratio, but as this aspect ratio is further docroasod, 
i^iCp drops sharply. In this caso, at aspect I’atios of approxi- 

mtely 0.5, becom.oa zero and assumes large negative valties 

with further reductions in aspect ratio. For applications to vei*y 
small aspect ratios, however, tho theory may require modificatiaas. 
The data in figure 3 show that as tho Mach number is increased, tho 
aspect ratio corresponding to zero valuo of becomes smaller. 

Tho tip effects shown in figure 3 fox* tho wing of 45° sweep- 
back are similar for other wings of different swoopback at appro- 
priate aspect ratios. The conversion formulas indicated in tho 
section entitled "Analysis" indicate that tho aspect-ratio effects 
for wings of difforent sweopback correspond qtialitatively for equal 
values of tho aspect-ratio paramotor A tan A* Tho Mach numbors 
for each of the wings differ, however. An aspect ratio of 0.8 
(fig. 3) for the wing of 45° sweepback at a Mach number of 1.10 
for example corresponds to an aspect ratio of 0.8 cot A at a Mach 

number equal to \jl + j(l.l0)^ - ij tan^ A for any other ■vlng of 
sweepback angle A • 

Variation of wing-drag coefficient- with Macli number, swoepback, 
and Blendenaes a ratio .- Figure 4 ahovs the variation of Cp with M' 
for different sweepback angles with constant slenderness ratios. 

The slenderness ratio represonts the ratio of the wing semispan 
measured along the leading edge to the maximum ttiicknesa of the 
center section. Tho data are presented for sweepback angles of 30°, 
45°, 52.5°, and 60° with slendemesa. ratios of 25 and 50. The 
wings in figure 4 for the difleront slondornoss ratios and s'vroep- 
back angles are assumed to have tho same ^d.ng area end tho salhe 
profile normal to tho >rfng loading edge. Tho slondcmcss ratios 
aro baaod on a thickness ratio of 0.10 measured in a direction 


r 
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normal to the "wing leading edge. Ihe thiclaaees ratio t/o measvired 
in the fll^t direction, therefore, varies "with sweephaci: as cos A 
or is equal to 0.1 cos A. The gispect ratio is reduced .vith sweep- 

hack hy the factor cos^A The aspect ratio is related to the 
slenderness ratio hy the formula: 

^ cos^A 

The plan forms for the different wings are shown in figure 4. 

The results in figure 4 show that, in general, the drag coef- 
ficient decreases with increasing sweephack* At Mach numhers for 
which the l-Iach lines are appreciably ahead of the wing leadi n g edge, 
increasing the slenderness ratio or aspect ratio gives inqportant 
reductions in calculated wave-drag coefficient. At Mach numbsps for 
tdiich the Mach linos approach the wing leading edge (M — j-sec A), • 
however, short wj.de wings give appreciable reductions in wave -drag 
coefficient. The figure also indicatos that 1die hi^ost wave-drag 
coefficients for tho normal range of aspect ratio occur at a Mach 
number equal to sec A 

Effect of aspect ratio on wing-drag coofflcient .- Figure 5 
Indicates the effect of aspect ratio on tho wavo-drag coefficient 
for the wing. The data in figure 5 show tho wavo-drag-coefflcient 
Cp tan A 

parameter — plotted against the aspect-ratio parameter 

100 (t/c)- 

A tan A These res ults a re shown for various valtios of the Mach 

dumber parameter - 1 cot A, which corrospond to a range <jf 

Mach number from 1.0 to the secant of the angle of sweepback, 
or 1 ^ M < soc A. 

Figure 5 shows that for a given value of the Mach number 
parameter, tho maalmuni wave-drag coofficlont occtcrs at a definite 

aspect ratio. For ezamplo, if A = 45° and - 1 cot A = O. 3 IO 

that is, M = 1 . 05 ) the maxiTmm value of Cjj occurs at an aspect 
ratio of 0.^. If tho aspect ratio is decreased to values smaller 
than 0 . 85 ^ Cp drops very sharply. Similarly, as tho aspect ratio 

is increased frota 0.8^, also decreases. Thizs, in general for 

the Mach number parameter corresponding to \^~ - 1 cot A =* O. 3 IO 
the znoTl . TmTm value of Cjj occtirs at an aspect ratio equal to 

0.^ cot A 

Application of curves of figure 5 bo wings of arbitrary B\ 7 eojp- 
baok and aspect ratio .- The scale labels and curves of figure 5 
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apply to the series of-vlngs that may he derlvecL frcsa a hasic 
lO-percent-thlck, 45° swept-hack ving. The lahols express the 
transformation equations (I3) to (15)* If A is sot equal to 45° 
and t/c is set equal to 0.10 in these labels to correspond to 
the basic wing, the ordinates become simply Cj), the abscissa A 

(aspect ratio), and the curve parameter ~ 1. The results in 

figure 5 inay be applied to all ewept-back wings covering a range 
of aspect ratio from 0 to 10 cot A corresponding to a range of 
Mach number from 1.0 to sec A* The data apply specifically to 
untapered irfjigs wltli biconvex pi'ofiles at zero lift with Idle wing 
tips cut off in the direction of fli^t. The I’esults, however, may 
be applied to indicate approximate results for profiles similar to 
the biconvex. 

The following example is given in order to illustrate the use 
of figure 5. I’or a given wing 


A = 70° 

A r, 3 - 

- a 0.08 
0 

M = 2.20 

In order to find. 0^; 

A tan A = 8.24 
\llf - 1 cot A = 0.715 

From figure 5, for A tan A = 8.24 and \/m^ ~ 1 cot A = O.715 

Ct) tan A 
-=0.0123 

100 


r 
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therefore 


C 33 = 0.C0286 

CONCinDUG Rfl^IARES 


A theoretlceG. investigation has "been made of the supersonic 
wave drag of untapered swept-back wings at zero lift. The wing 
sections were biconvex and the wing tips were considered to be cut 
off in the direction of fli^t. The investigation was limited to a 
range of stream Mach nmber from 1.0 to a value corresponding to 
that at which the Mach line coincided with the wing leading edge . 
For tills range of Mach number, the following conclusions have been 
drawn: 

1. In general, the calculated wave-drag coefficient decreased 
with Increasing sweepback. 

2, At Mach numbers for which the Mach lines are appreciably 
ahead of the wing leading edge, increasing the slenderness ratio 
or aspect ratio gave inmortant reduction in the calculated wave- 
drag coefficient. 

3» At Mach numbers for "vdiioh the Mach lines approach the wing 
leading edge (Mach numbers approaching a value equal to "the secant 
of the angle of sweepback) , decreasing the slenderness ratio or 
aspect ratio reduced the calculated wave-drag coefficient. 

4. The highest calculated wave-drag coefficients for the 
normal range of aspect ratio occurred at a Mach number equal to 
the secant of the angle of sweepback, 

5 . The maximum wave-drag coefficient occurred, at a definite 
aspect ratio which is determined by the Mach number and angle of 
sweepback. 


6 . For aspect ratios greater "than 1 /vM^ - 1, where M is 
the Mach number, the increment in wave-drag coefficient for ihe 
wing contributed by the tip was zero. 

7 • 5he variation of the drag with Mach number obtained for 
one sweepback an^e for appropriate aspect ratios may be presented 
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in a Tinifled form so that the drag for the complete range of sweep- 
haclc engle^ aspect ratio, and. Mach number may be directly determined, 
from a single chart. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December 17, 19^6 
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APPEiroil A 

FORMULAS FCE EIE ZNT23GRA1TD u EQUATION (3) ‘FOR FUflTE 
UHTAPERED SWEEP -BACK WINGS OF LICQN^'EX ERCFILE 
AI ZERO LIFE 


m = cot A < 


= P 


In. order to aatiafy the houndary condition for a finite swept- 
hack wing of biconvex profile, the integrand u in equation (3) 
may he expressed in terms of components caused hy elementary source 
linos as follows: 


u = 


^ c ,0 ^,0 ^ c ,0 


" (^,0 “ 5 ^ c ,0 + |« 0,0 " 5^,0 


^/m,h ~ -h^ 


(Al) 


where the subscript notation indicates the origin of the source 
line. The bars over u refer to the source lines caused hy the 
opposite wing panelj that is, u indicates a source line with a 
i-evorsal in the sign of m. 

In equation (Al), tlae u-functions are given by the real parts 
of the following expressions: 


u (x. 


y) 


= I cosh 


■1 


X - g - m^^(y - T]) 
p|y - ^ - m(x - g) 


u (x, 


j) 


-1 ^ + mp^Cy 

= I cosh 


n) 


P |y - 7} - 1 - m(x - g) j 
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where t| repreBents the origin, of the olementary source linos « 
For the "biconvex profile. 


dx c 


and frcau equation ( 5 ), the f so tor 


2t V 
c It 


m 






The symbol i in equation (Al) refers to an integration opera- 
tion which represents the influence of a uniform distribution of source 
lines along the chord of the biconvex profile beginning at the 
position 2,^ T|, Ohls symbol is defined by the following espiressl one : 





Ix-pjy-Ti 


^ cosh’i 

dg' pjy - T) - m(x - ^')| 


dS' 


= (y - n) 



i Ni 1 - nrE^ 

dx 

L “ 


cosh'd 


X - 5 

3!y - n 



m(x - g,) 

y “ n 


cosh”^ 



- ^ 

^ m(x -^i) 


y - TJ 

mmJ 


Ca2) 


whore £' is iho variable of integration representing the ^ -coordi- 
nate of the origin of each source line in the distribution of source 
linos . 


r 
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For a liioonTex profile 


_ _ Y m d^z 



Y m 
^ ^1 - 


Eq.'uation. (A2) may 136 eaqjresaed as a fxuictlon of 


X ^ 

that 1b, 

y - T) 


7) . - n) ^(Trf ) 


Then 


y) =-{y - f,) f 




- (y - Ti) 


Tho limitH of IntsfpiratiOns vith regard to x for the section 
drag coefficients and with rfegard to y for tho total wing-drag 
coefficients are dlscnssed* The u-dtmponents catised "by each of the 
elementary source lines are aero at all points outside of tho 
respective Mach cones. The functions for the u-intogrand in eq.ua- 
tion ( 3 ) are therefore evaluated along the section for values of x 
"beginning at the forward "boundary of the Mach cone. This integra- 
tion gives tho section-drag-coefficlent cemponents. Similarly^ in 
order to obtain tho wing-drag coefficient, the section drag coeffi- 
cient components obtained from the respoctive u-functions are 
evaluated along the \Jing span for values of y contained within 
the Mach' cone. The followirig table refers to one side of the wing 
(x and y positlvo) and shows the limits of integration for x 
and y for the roq.tiired u-functions: 



Hlfi ! pi 
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Limits of integration 


u-coccponents 

X 

7 

Lo-wer 

limit 

Upper 

limit 

LoTOr 

limit 

Upper 

limit 

'\),0 ^,0 
1 

y/m 

1+ c 
m 

1 

0 

i 

i 

h 

""c,0 

1 1 - 

5^^0,0 

i 

1 

y(3 + c 

i 

Z + c 

m 

. . 

0 

1 

h 

1 

^h/m, h 
1 

|(mp + 1) - yP 

1 

! 

^4 C ' 

m 

j 

1 

1 

, mo 

h •* 

mp + 1 

(if A > ^ 

\ mp + 1 / 

0 flf A < — 

y “ Dip + ly 

li 

1 1 
Af. 

1 

yP + ^(m3 + 1) 

1 : 

1 

1 

^ 4 C 

m 

i 

! 

1 

* • 1 

li(rap + 1) mo 
1 - mP 

(±f A > — ia— \ 

\ mp + 1 y 

0 (if P.< ^ \ 

\ ='mp + 1 / 

• ".— 1 

h 


I 
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APEEUDIX B 


FOEMULAS FOR WAVB-niAj3- COEFFICIENTS FOE FINITS UNTAFEEED 
SWEPT "BAaK: WESGS OF SYI>fi'IETRICAL BICOJJVEX EROFUE 
AT ZERO LlbT [m = cot A < i 

L ^ 3 J 

Sectioa Drag Coeff iciontc 


In the following analyaie the quantities y and K are 
employed nondimensionelly in terms of the seinl chord. The eq.uations 
for the drag coefficients in all oases refer to the real parts of 
tlie indicated expressions . 

Section drag coefficients wiidiout tip effect .- The section drag 
coefficient for the given wing at a spanwise position y and Mach 
number M without the tip effect was found to be as follows; 


The term 


2m. 


represents a convenient integration liinit 


1 - mp 

which indicates the intersection of the Mach line from the center- 
section trailing edge with the wing leading edge. For y = Em < ^ 


(y) = “ (I) fooBh'^ _ 2 coBh“^ 

“•CD It \c / \^ \ 


3\/l - m’^ I 


2m' 


- a(aK3 - as - 1) coah-^ * g 

2m‘(E + 1) 


+ 4k(2K“ - 3 ) cosh ^ ^ - ■ -}- K®[o!:(l - m'^) 

- \jil - m’^) [(K + 2)2 - (Kk’)2j 


1 - mp 


(Bla) 
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E’er y = Km > — ^ the following expression is addncL to eq.ua- 
1 - m3 

tion (Bla) . 


8 /t\^ „ ,„3 


rt V. c/ 


. -1 E - 2 
m <,E cosh 

I 


3^1 - m'^L 


2 ooeh 


-1 - m*^) - 2 


2m' 


^ - N , -1 k(i + m j - 2 

+ 2 v^ciE^ - 3h + 1 / cosh ■ 

2m’(E - 1) 


1 


Kr^(l - m'2).[(K - 2)2 - (Km')2] j | 


(Bib) 


where m' = mp. 

Foi’ the special case m a — — . the Mach line coincides 

\/m2 - 1 

with the wing leading edge, and the expression for the drag coeffi- 
cient obtained as a limiting case (m' — j- 1) for all values of y 
becemes; 


w (i) 

» IC \ 0/ 


'd 


. i3 , -1 y ‘ + 2 2 
m <y' cosh * 

y ' 3 


y'(3y' +-4)(y' - 1 ) - 2 


\/y' + 1 


(B2) 


•rfiere 


y' “ yP 


At the center section, ■vdiere y or E = 0, equation (Bla) bocamesi 
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,.-1 1 

■ ■■■ - cosh, — 

L 2^ mp 

/I - m p 


At the center section, for y = 0 and m = i, eq.uation (B2) 
■beccsmos : 


- ^ /if 

Vc/ 


IncreBient In section dra^ coefficient caused b;r ving tips • “ 

The Increment in c^ catised hy the tips depends cn various factors, 

such as the sweep angle, eapeot ratio, and Idach nimiber. The following 
types occur in an untaperod vlng: 

If aspect ratio of the wing is eq.ual to or greater 

than i/\/m 2 - 1 each tip affects solely its own haiJT of the wing. 

In this case the effect of the tip is limited to tho region of the 
w5.ng outboard and roar5^ard of the front Mach liru 3 oa-iginatlng from 
this tip. (See fig. 2(h).) The region of the wing affected is 


between values of 


from h ^ 


mP + 1 


to h. 


The increment in section drag coefficient at a Mach number M 
and spanwise position y caused by the tip was found to be aa 
follows ; 




26 


M-ACA TN No . 1319 


where the subscript a iuMcates that the x-axis is shifted to the 
tip section, and jg.' = and m' = In the plan form of the 

wing 


y = h + y. 


In equation (B3) values for y^ may be taken from - 


2m 


m(5 + 1 


to 0. When the Mach line is coincident iri.th tlie leading edge of the 

the expression for ^®cames: 


1 


L • 

\%> 

■, 1 

t f 

^a 

— ) m 
c/ 

- 

\jya‘ + ^ 






cosh 


•1 ya' + 2 


y® 


3 . 


(3Ta'^ - 25-a' + **) 


n. If the aspect ratio of the wing is less than 


(BJO 

the 


Vm2 - 1 

tip on the opposite wing contributes an increment in c^ in addi- 
tion to that discussed under type I. The increment in at a 

section caused by the opposite tip was obtained Jja the following form: 



.2 


~ — (l4+ m''^) - 3h'y, ' + 2h'2 
12 ^ 


cosh 


„1 y-b' - 2(h' -m') 


m'y^« 


- m'2 

-1 y-b ' t 


cosh 


2m 



+ m 


’^-6y^‘2j^' + 3yb'(2h’2- 

- 2(h' - m») 

■ - h’ + m'^ 





(B5) 
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Tdaere the eubscjript h Indicates that the x-aacLs is shlfijed to the 
opposite tip section, and -vdiere m’ = mp, and h* = hp. 

!Dn. the plan farm of the wing 


7 = 7^ 


- h. 


The limits for y-u to he used in equation (B5) depend on the 
value of aspect ratio A relative to the parameter Thtis 


mp + 1 

(a) If the aspect ratio of the wing is greater than 


2m 


mP + 1 

the front Mach line frcm the opposite tip intersects the trailing 

edge at a value of y-„ = so that values for y, in 

“ 1 - mp - ^ 

Ofy, 

equation (B5) may be t^en from ~ ^ to 2h at the tip. 

. 1 - mp 

(See fig. 2(b).) 


(b) If the aspect ratio A is equal to or less than 


2m 


mp + 1 

the front Mach line frcm the opposite tip intersects the center 
section and values for y^ in equation (B5) may be taken frcm h 

to 2h. In this case, the increment in Acg^ di.sctsesed under 
tjrpe I is obtained at spanwise positions of y^ from -h to 0. 


When m = 


- 


:, equation (B5) becames: 


(if ^ ’ 2)\/7^’^2(h‘-l)2-y^'2 

-(|7t’^’3h’y^' + 2h'2-i^ ^osh~^ 

^ 2 [^•^“^7b*^* + 3y^*(2h«^-l)~2h»^ + 3h«-l]^/r^1 
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Tho total increnoent In drag coefficient at a section caused hy 
the tips is given hy 


^d = 


and the total drag coefficient at the section is 


°d = ^'d^+ 


Wing-Drag Coefficients 

Wing-drag coefficient Tdthout tip effects.- If the aspect ratio 

pjj. ’ 

is equal to or less than , the wing-drag coefficient without 

1 - mp 

the tip effect is 


. 8/tf^ A'2 /, .-lA‘ + 2m» 

[l2m'3\ 


A’m‘ 


- 6a' cosh 


^ ~ - \|A'^(i - m*^) + hm’(A' + El*) 
m’ 


+ 2A’\jl - )+ i - 


co.h-1 

2m’^ 


(an '3 t 3 A'm'= - a' 3 ) cosh-^ 


2m*(A’ + m') 


+ (aA*^ - 6A’m'^) 


cosh 


-1 1 + m 


i2 


2m' 


where A’ =» A3 and m* => m3 
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If the aspect ratio is greater than 
coefficient -without tip effect is. 


2m. 


.1 - 


-bhe -wing-arag 


Cj) = ■§ (it m 
•^00 rt \c/ 


.12m 


,3 


/ ^.-lA•+^m* .-lA'-2m‘ 

3A* (cosh + cosh 


A'm' 


A'm* 


- 6A' cosh" 


1 i_ - uA'^fl “ m'^) + l«a’(A' + m‘) 
m' ' 

- \|a'^(i - m'^) + - A') + SA'^ - m’® 

.3 A-(l - m’^) + 2m’ 


3m‘ 


3^1 . aji 2 


2m 


,2 


- cosh 


-1 A*(l - m'^) - 2m' 


2m 


.2 


* (sA'm'S - 2m'3 - A'3) coeh-^ 

2m’(A* - m') 


i2 o_i3 « ,3'\ , -1 A‘(l + m'^) + 2m’ 

+ \^3A'm' + 2m’ - A’~’/ co£dx — i 

• 2m’(A’ + m') 


+ (2A'^ - 6A’m*^) cosh"^ 

• 2m' 


Increment in -wing-arag ooeffloient caused >y -wing tips «~ If 
the aspect ratiQ of -the wing is equal to «r greater than l/p, the 
contrittEbion of the wing tips to -fche -wing-irag coefficient is Jtaro. 
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Xf tlie aapect tatio of the ving is equal to or greater than — , 

-hi 

the total Increscent in Cp caused by the local tip, or is 

seroj and the incroiaeBt Z£!p is obtained solely Crcaii the effect of 
the tip located on the opposite half of the ving, or 
this case, integration of equation (B5) over the uing yields: 




L3m»3jl - a '2 


(A’ - aii')(A’ + m’)^ cosh'^ . t . i: 

A' + m* 


+ (a’ + 2m') (A’ - m')^ cosh 


3m' 


■1 1- n'A' 

I*' - m'l. 


L»2 _ ^i2^2 + m'^) 

“1 

cosh”^ i- 


1 - A-" 

L ^ \ 

J A' 

3m'V 

J 


If the aspect ratio of the ving is less than 


2m 


the 


mP + 1 

increment dCp is affected by both wing tips. For this case, 
integration of equations (B3) and (B5) yields: 


^ - ! (If 


m 


.31a' 






(A* - 2m') ( a' + m’)^ cosh 


,-l m 'A ' + 1 
A' + m' 


+ (a' + 2m*)(A' - m')^ oosh"^ , - 2m'^ cosh ^ 
- (sm'^ “ 3m’^A’ + A'^) cosh"^ .'.A'f?- •*• 


|2m’(m* - A')( 


+ ~ I 6m’2 -A'^(2 + a'2)^ cesh’^ - A') + A‘^“(l - m'^)j 


3m? 
A 


^ ooBh-i - al,£Ti:2'\ 

4m'3 “*A‘ 3m'^ J 
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Total •wlng-'drag coefflclen.t »- The total "wing drag is obtained 
as the sum.. 


9d * ^ 

CO 


vhere the ccanpanenta said aa?o oalculated from the fore- 

going equations for lihe ning-drag coefficiont approixriate to the 
aspect ratio of the ving. 


Wing-drag coefficient for special case m = ■=»- When the Mach 
& 


line coincides vith the -cdne leading edge 


(” - 1 > 


the -wing-drag 


coefficient obtained for all aspect ratios as a limiting case (m* 
is equal to the real part of -fche folio-wing expressicn: 


1 ) 



Wing-Drag Gooff icients at Mach Number of 1.0 

♦ 

The drag coefficient for the -wing at M = 1.0 may be expressed 
in teims of the f oUo-wlng forraulaB -thich wore ob-balned by integrating 
over -fche -wing -bhe limiting values for in equations (Bl), (B3), 

and (B5) as the factor p = approaches aero: 
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than 

Cjj =r 


(1) If the aspect ratio of the wing is eq.ual to cr greater 


Sm 


L 

3« 


(if sM “ [(-K^ SE + 4) log, (K 


+ 2 ) - 3K ^ log ^ 2K 


+ (kK^ - 6K - 2 ) log (K + 1) + K^(K - 2)] 


as: 


ph/m 


m 


1)2 


[(-s:^ + 6 E - 4) log^ (K - 2 ) - (k^ - - 4) 108q(K + 2 ) 


+ (4e: 3 > gjj.) log (k^ - 1 ) + 2 log SjLi - ,^3 ^ 

B ® E + 1 ® J 


as 


'h 

(jh-2m LL 


(y “ h)3 3 (h - y) _ p 

2m3 ^ 


l^log^ 2 (y - h + 2m)J 


- j - 


(y + h)3 ^ 3y(y^ - h^) _ 3 (h + y) 


2in? 




m 


m 


Jlogg m(y + h) j 


+ -(7y - 3h - 2m)(y - h + 2ia) ^ dy] 

4m '^ 


(B6) 


y 

m 


■where K 
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In equation (b 6), the first tw> integrals represent the arag 
coefficient for the ving witho^it the tip effect; T&ereas, the last 
integral represents the effect of the tips Far this case, 

in Tdiich the aspect ratio is equal to or greater than Sm, the tip 
effect ACj)^ is zero; hence tJie integral for the section increments 

in c^ is not given in equation (B6) . Equation (Bo) lias heen 

solved for a sweephack angle of lj-5° and the results for this sweep 
angle may he converted to other w?.ngs of arhita’ary s'vreephack hy the 
formulas ( 13 ) to ( 15 )* ^'or the wing of ij-5° sweephaok, ‘m = 1 
and' A >2, equation (b 6) yields the following result: 




(-A^ + 12A + 16 ) logg (A + 2) 

(-A^ + 12A - 16 ) log^ (A - 2) + 2 a(a^ - I 2 ) log^ A ~ Ii-aJ 


(B7) 


(2) If the aspect ratio of the wing is amalleor than 2m, the 
upper limit of the first integral in equation (b 6) is reduced from 2 
to h/m, the second integral vanishes, and the lower limit for the 
third integral is reduced from h - 2m to zerot Par this case, 
however, in which the aspect ratio is less than 2m, is not 

zero, and the following integral must he added to those in equa- 
tion (b 6) to obtain C^: 



Tdiere 


y^ = y - h 
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Por the vlng of 45® sveophaclc, m a 1, and A < 2, ectuationa (b 6) 
and (b 8) yield the sanie restilt for Cp ajs that obtained for values 

of A greater than 2, aa ezpreesed by eq.mtlcaa (B7) • In this case, 
the reeil part of logQ(A - 2) is used. 
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Section rtrag coePP/c/ent ^ 


Fig. 2a 
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Mach lines 
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Distance Prom center, y, scmichords 
ca)Wiihoui ftp effect a! several thach numbers. 
Figure 2.. ^Tgpical distnbcitioneicf section wgye-drai 


difectiorFo.lO. 
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Fig. 



(b)'^jfh ftp effect. Aspect rat to j f.SQj Mach numberj I.IQ. 


Figure 2.. -Concluded 
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Aspect ratio ^ A national advisory 

CONHITTEE FOR AERONAUTICS 

F tgareS.— Typical va nations with aspect ratio oF the mcrer- 
/nents /n win^ wove-dra^ cOBfficienf caused £>y '*Jmq 
tips at d iff&rent tAach numbers^ Biconvex pi^of//& at 
xero hft I svieepbact angle j4S°j thickness ratio m flight 
direct lonj O-tO . 






Drag-coePFiQieni pararneier 
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